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We analyzed localized charge time evolution in the system of two interacting quantum dots (QD) 
(artificial molecule) coupled with the continuous spectrum states. We demonstrated that Coulomb 
interaction modifies relaxation rates and is responsible for non-monotonic time evolution of the lo- 
calized charge. We suggested new mechanism of this non-monotonic charge time evolution connected 
with charge redistribution between different relaxation channels in each QD. 
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I. INTRODUCTION 

QDs are unique engineered small conductive regions 
in the semiconductor with a variable number of strongly 
interacting electrons which occupy well-defined discrete 
quantum states, for this reason they are referred to as 
"artificial" atoms [HJl]- Several coupled QDs form an 
"artificial" molecule [Ij],[H and can be applied for elec- 
tronic devices creation dealing with quantum kinetics 
of individual localized states [H , [fil , Q . That's why the 
behavior of coupled QDs systems in different configura- 
tions is under careful experimental [8[ , [9j and theoretical 
investigation (Toj.(Tl|. It was demonstrated experimen- 
tally that coupled QDs can vary from the weak tunneling 
regime (coupling with the leads is smaller than interac- 
tion between the QDs) to the strong tunneling regime 
(interaction with the leads exceeds the QDs coupling) 
3], [12]. One of the most perspective technological goals 
of QDs integration in a little quantum circuits deals with 
careful analysis of non-equilibrium charge distribution, 
relaxation processes and non-stationary effects influence 
on the electron transport through the system of QDs 
[HI > [3 > [HI > [3 , [13] • Electron transport in such systems 
is governed by Coulomb interaction between localized 
electrons and of course by the ratio between the tunneling 
transfer amplitudes and the QDs coupling. Correct inter- 
pretation of quantum effects in nanoscale systems gives 
an opportunity to create high speed electronic and logic 
devices [H| , [l9| . So the problem of charge relaxation due 
to the tunneling processes between QDs coupled with the 
continuous spectrum states in the presence or absence of 
Coulomb interaction is really vital. Time evolution of 
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charge states in a semiconductor double quantum well 
in the presence of Coulomb interaction was experimen- 
tally investigated in [20(. Authors manipulated the lo- 
calized charge by the initial pulses and observed pulse- 
induced tunneling electrons oscillations which were fit- 
ted well by an exponential decay of the cosine function 
and a linearly decreasing term. Time dependence of the 
accumulated charge and the tunneling current through 
the single and coupled quantum wells in the absence of 
the Coulomb interaction were theoretically analyzed in 
21 1, [22j], [23[. But the authors took into account only 
two time scales which determine charge relaxation and 
neglected the third time scale which is responsible for 
charge redistribution between different quantum wells. 

In this paper we consider charge relaxation in a single 
QD and double QDs due to the coupling with the con- 
tinuous spectrum states. In the case of two coupled QDs 
tunneling to the continuum is possible only from one of 
the QDs. We have found that on-site Coulomb repulsion 
even in one of the dots results in significant changing 
of the localized charge relaxation and leads to formation 
of several time ranges with strongly different values of 
the relaxation rates. We pointed out that the leading 
mechanism of non-monotonic charge relaxation is charge 
redistribution between the relaxation channels in one of 
the QDs due to the Coulomb interaction. 



II. NON-STATIONARY TUNNELING 
PROCESSES IN THE SINGLE QD 

First of all let us consider QD coupled to an electronic 
reservoir (conduction electrons states have energies £&). 
We assume that the single particle level spacing in the 
dot is larger than all other energy scales, so that only one, 
spin-degenerate level of the QD spectrum is accessible e\. 
Such a system can be described by the Hamiltonian: 
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FIG. 1: Solid line: Localized state filling numbers time evolution in the QD with energy level ei when distribution function of 
the continuous spectrum electrons is taken into account. Dashed line: Filling numbers time evolution in the QD with energy 
level £i without contribution from the conduction electrons many-particle effects. Dash-dotted line: Filling numbers evolution 
only due to the many particle effects caused by the presence of the distribution function. Tunneling transfer rate 71 = 0.3 has 
the same value for all the figures, a) ei = 1.3, b) ei = 0.3, c) Ei = —1.3. 
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where T k is a tunneling amplitude between the QD and 
the continuous spectrum states which we assume to be 
independent of momentum and spin. c\jc\ and c^/cv 
electron creation/annihilation operators in the QD lo- 
calized state and in the continuous spectrum states (k) 
correspondingly. 

Let us assume that at the initial moment all charge 
density in the system is localized in the QD and has the 
value njfO) = uq. We shall use Keldysh diagram tech- 
nique [241 ] to describe charge density relaxation processes 
in the considered system. Time evolution of the electron 
density in the QD is determined by the Keldysh Green 
function GJ\ which is connected with the localized state 
filling numbers in the following way: 



G&(t,t) = tni(t) 



(2) 



System of integro-differential equations for the Green 
function G^,(t, f ) has the form: 



where continuous spectrum states Green function 
GW(t, t) and inverse localized state Green function 
G^i _1 m the absence of tunneling processes have the 
form: 



G° k %(t,t) = -i<d{t - t )e- iek{t - 1 ^ 
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In equations © and Q integration over intermedi- 
ate time arguments is performed. Finally the solution of 
equation ((4]) can be written as: 

G<{t,t) = m(o)e-*»* + y; f [ e(t-ti)G(t-t 2 )- 

k J° ^° 

■ dtidt 2 f(e k )e~ lEkitl - t2) . e -«(*-*i) e ^~i*(*-*2) 

(6) 

where we define 
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and consequently one can obtain the following equation 



^-continuous spectrum density of states which is not 
a function of energy, /(w)-Fermi distribution function. 

Performing integration in expression (|6]) and replacing 
summation over k by integration over u) one can get final 
expression which describe filling numbers evolution in the 
quantum dot due to the interaction with the continuous 
spectrum states: 
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In general, localized charge relaxation law differs 
from the simple exponential law even in the absence of 
Coulomb interaction. Similar expression was obtained 
for the initially empty localized states time evolution by 



means of Heisenberg equations in [21[ . 

Figure [1] demonstrates the localized state filling num- 
bers n\ {t) time evolution for the different initial positions 
of the energy level in the QD. When the continuous spec- 
trum electrons have Fermi distribution function, charge 
density relaxation law strongly differs from the exponen- 
tial law, especially when condition \e\ — £f\ < 7i is valid 
(solid line in Fig|TJ). This difference can be seen even 
when t < r^-j^q if condition |ei — ep\ ;g> 7! occurs. 
It is clearly evident that when contribution from many- 
particle effects in the continuous spectrum states is ne- 
glected charge relaxation demonstrates simple exponen- 
tial law (dashed line in FigfT]). Contribution only from 
the continuous spectrum many-particle effects is depicted 
by the dash-dotted line in FigfTJ 

Stationary distribution can be achieved for t — >• oo: 
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FIG. 2: Schematic diagram of energy levels in the system of 
two coupled QDs. Second QD is also connected with contin- 
uous spectrum states. 



G?f(i-i' 



-ie{t-t')e- iei 

-iS{t - t ')e-**(*-*')-r(*-*') (12) 



where 7 = KVuTu is the tunneling relaxation rate from 
the second QD to the continuous spectrum states. 

Retarded electron Green's function Gf x yields density 
of states in the first QD and can be found exactly from 
the integral equation: 



III. NON-STATIONARY TUNNELING 
PROCESSES IN THE SYSTEM OF COUPLED 
QDS 

Let us now investigate charge relaxation processes in 
the system of two coupled QDs with single-electron en- 
ergy levels £1 and £2 correspondingly (FigfS]). QD with 
energy level £2 is also connected with the continuous 
spectrum states. Hamiltonian of the system under in- 
vestigation has the form: 



QRrr2 n 0Rr<R 



(13) 



The eigenfrequencies E ly2 of equation (fT3"|) are deter- 
mined in the following way: 



(E-£ 1 )(E-£ 2 +n)-T 2 = 



Ei, 2 = -(£1 + £2 - »7) ± 2^( £ i -£2 + i7) 2 + 4T 2 (14) 
Finally retarded Green's function can be written as: 



a a k,cr 

+ J2 T ( c t* c 2<? +4x c i^) +J2T k (c+ a c 2a + 4 <J c ka )(ll) 

(7 k,(T 

T and T k are tunneling amplitudes between the QDs 
and between the second dot and the continuous spec- 
trum states correspondingly which we assume to be in- 
dependent of momentum and spin. C\ f c\{c\ j cq) and 
c^T/cfc- electrons creation/ annihilation operators in the 
first (second) QD localized state and in the continuous 
spectrum states (k) correspondingly. 

We assume that at the initial moment all charge den- 
sity in the system is localized in the first QD and has the 
value ni(0). First of all we have to calculate exact re- 
tarded Green functions of the system. In the absence of 
tunneling between the QDs Green functions G\^{t — t ) 
and G^it — t ) are determined by the expressions: 



G*(t,t) = -i&(t-t)( 



E 1 -e 2 +i'y c _ El(t _ t ' ) 
Ei — Ei 
E2 — £2 
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Let us now analyze time evolution of the electron den- 
sity in the considered system. Electron density time 
evolution is governed by the Keldysh Green function 

G<(t,t) 



Gf 1 {t,t) = in 1 {t) 
Equation for Green function Gfi has the form: 



(16) 
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Acting with operator G°f 1 it can be re-written as: 



it has the following form: 



G™- 1 G<(t,t') = (i?- t -s 1 )G<(t,t) = 

= T 2 dt 1 G *(t,t 1 )G< 1 (t 1 ,t) + 
Jo 

poo 

+T 2 / dtiG^it^G&fat) 
Jo 

or in a compact form: 



(18) 



Gfi(t ,t) =iAe^ Elt+lE ^' + iB e - iElt+iE ^' + 

+ iB*e- lE2t+lE '^' +iCe- lE2t+lE t t ' (24) 

As far as solution has to satisfy homogeneous integro- 
differential equation (TT91 (without right hand part), after 
substituting expression (|22|) to equation (p~9]> one can find 
the following relation: 
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Green function G^ityt) is determined by the sum of 
homogeneous and inhomogeneous solutions. Inhomoge- 
neous solution of the equation can be written in the fol- 
lowing way: 



G<(t,t')=T 2 f dti / dhG^it-h)- 
Jo Jo 

"^22 (*1 — h)G^ 1 (t2 - t ) 



(20) 



Green function G® 2 can be found from equations ^ 
and ((H) substituting index 2 instead of 1. 
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Using the initial condition: 

G<(0,0) = m? 



(25) 
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Time dependence of the filling number n\ (t) in the first 
QD can be obtained: 

m(t) = T»S • (A'e-*^ 1 - 8 ^* + 2Re(B'e- t(El - E ^ t ) + 
+ d e-^*-^) (27) 

where coefficients A', B 1 and C" are equal to: 
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(21) 



Green function G®^ contains part with exponential de- 
cay, oscillating term and part which determine station- 
ary solution. In what follows we'll consider the situation 
with £i ~ e - F >> 1 for simplicity. It means that stationary 
occupation number in the second QD in the absence of 
coupling between QDs is of the order of — - — << 1. So 
we can omit corresponding terms in expression (|21[) for 
function G^ ■ 

As we consider initial charge to be localized in the first 
QD, Green function G2^(0, 0) ~ and Green function 
Gfi (t, t ) can be determined by the solution of homoge- 
neous equation. Homogeneous solution of the differential 
equation has the form: 



G<(t,t) = h(t)e- iE ^ + Ht)e- iE * (22) 
Since G < (t,t ) satisfies the symmetry relation: 



(G<(t,t))* = -G<(t,t), 



(23) 
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(28) 



Electron density time evolution in the second QD is 
determined by the Green function G 2 < 2 (i,t ) with initial 
condition G^O, 0) = 0. Green function G^itit ) can 
be found from equation similar to equation (|19[) with the 
following indexes changing (1 ■<-» 2). Due to the initial 
conditions 712(0) = 0, rti(0) = no, filling numbers in the 
second QD n 2 {t) are determined by the inhomogeneous 
part of the solution. Electron filling numbers time de- 
pendence in the second QD n 2 (t) can be written as: 

n 2 (t) = (£> e -i(£i-Er)t + 2Re(Ee- liEl - E ^ t ) + 

+ Fe-'^-m*) (29) 

where coefficients D, E and F are : 



D = F = 



T 



\E 2 -E l \ 



E=-- 



T 2 



\Eo 



(30) 



Expressions (|27|) and (|29[) looks like there are three 
relaxation channels with different time scales. The first 
and the second relaxation channels are connected with 
relaxation rates (\E\ — E*\) and (\E 2 — E% \) . One more 
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time scale is connected with the expression {\E\ — E%\). 
This time scale is responsible for charge density oscilla- 
tions in the both QDs, when the following ratio between 
T and 7 is valid: T/<y > 1/2. 

In the resonance E\ ~ e 2 one can find four different 
regimes of the system behavior: 

1) Realization of the condition 2T < 7 leads to the 
absence of oscillations in the QDs charge density time 
evolution. In this case the following expressions are valid: 



7res — "fnonres — ^fres t \o V"^/ 

7 (£i-£2r 

It is not surprising of course that 7 res 3> Jnonres 

Let us take into account Coulomb interaction between 

localized electrons in the QDs. In this case interaction 

Hamiltonian can be written as: 



Ei 
Ex 



E* = 
E* = 



-i 7 (i - y /i - wjm 

-i 7 (l + y/1 - (4T 2 )/ 7 2 ) 

-vy 



2) When condition 2T <C 7 is fulfilled time evolution 
of the electron density in the first QD can be described 
by the expression: 



m(t) = n? 



2T< 



1 + 1 ' 



(31) 



In this case the main part of the charge decreases with 
the relaxation rate 



Hint = U 2 n 2a n 2 ^ a + U\ni a Tl\. 



(37) 



We shall use self-consistent mean field approximation. 
It means that in the derived expressions for the fill- 
ing numbers time evolution it is necessary to substi- 
tute energy level value £{ (i — 1,2) by the expression 
Si = Si + U- < rii{t) >. Then one should solve self- 
consistent system of equations for n,i(t). We shall analyze 
only paramagnetic case: = ni- a = rii. 

Such approximation can be applied when the following 
relations are fulfilled: 



\Ex~El\ «mm(|i;i|,|£; 2 |) 
\E 2 -El\^imin{\E x \,\E 2 \) 
lEx-EZl^minilExl,^) 



(38) 



lres = 2T 2 h 



(32) 



3) A special regime exists in the system when condition 
2T = 7 is valid. Relaxation of the charge in the QDs is 
non exponential: 



7 i)e 



-7* 



n 2 (t) 



(33) 



4) In the case when condition 2T > 7 takes place 
charge density oscillations can be seen in the both QDs 
with the typical frequency il = y/AT 2 — 7 2 , for 2T ^> 7: 



m(t) = n?e~ 7 * - [1 + cos(2Ti)] 



(34) 



Let's now analyze non-resonance case. If we are far 
from the resonance, relation \e\ — e 2 \ 3> 7, T takes place, 
and the filling numbers relaxation law in the first QD has 
the form: 



m(t) = n v 1 [ 1- 



2T 2 



(ei-e 2 ) 2 
2T 2 



g (ei-E2) 



T7* 



£2 



cos[(ei - e 2 )t] e-*] (35) 



Relaxation rates "f res and 7„ nres in the resonant and 
non-resonant cases are connected with each other by the 
relations: 



Inequalities (|38p mean that one can uncouple rapidly 
oscillating Green functions (G{\ and G^) and slowly 
changing functions ri\ (t) and n 2 (t) . Suggested conditions 
are analogous to the approximations which are used in 
the adiabatic approach. In the mean-field approximation 
the main effect of Coulomb interaction deals with the de- 
tuning changing between energy levels £\ and e 2 . So for 
simplicity we shall take into account Coulomb interaction 
only in the second QD coupled with the continuous spec- 
trum states. Further we'll demonstrate that the presence 
of on-site Coulomb repulsion in both QDs slightly modi- 
fies the obtained results. 

Let us discuss the application possibility of the mean- 
field approximation in the considered non-stationary 
case. For the stationary case when electron filling num- 
bers are changed by variation of the applied bias or gate 
voltage in mixed- valence regime [25[ conditions C//7 > 1 
and (si — £f)/i > 1 are important for the validity of 
the mean-field approximation. We are interested in the 
non-stationary effects, so these conditions are not so cru- 
cial. Slow variations of energy levels due to the localized 
charge time evolution allows to obtain reasonable results 
in the self-consistent mean-field approximation even in 
the presence of strong Coulomb interaction. We also 
want to point out that initial energy levels are situated 
well above the Fermi level (si — £f)/7 >> 1, so station- 
ary Kondo effect can not appear. Energy level in the QD 
connected with continuous spectrum is nearly empty. Of 
course non-zero electron density appears in the second 
QD during localized charge relaxation. So one can try to 
investigate non-stationary Kondo effect and estimate the 
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time scale of many-particle correlated state formation. 
The simple estimation of this time scale r is connected 
with the inverse width of Kondo peak. 



Kondo 



£2 



(39) 



Consequently relative values of the system parameters 
U/s2, U/j, T/7 and £2/7 determine it's behaviour. Our 
investigations deals with the typical parameters values 
demonstrated below on FigJSjFigJSl For these values of 
parameters jKondo ~ 1 X 1CT 2 7, so t ~ 1 x 10 2 • to, 
where To is the localized charge relaxation time in the 
second QD due to the interaction with continuous spec- 
trum. More careful r estimation can be based on the 
approach suggested in (2(| . We can consider energy level 
£2 position changing due to the effect of Coulomb inter- 
action U < ri2 (t) > to be similar to the influence of time 
dependent gate voltage on the electron energies in QDs 
leading to the spin flip. So following the logic of Glazman 
et. al. [2a| one can obtain: 
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FIG. 3: Schematic diagram of energy levels position for differ- 
ent signs of detuning in the presence of Coulomb interaction. 
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For typical values of system parameters (FigJSjFigin]) 

r -i „ icr 2 7 . 

Thus characteristic time of Kondo peak formation is 
much larger than the localized charge relaxation time. 



FIG. 4: Energy levels detuning time evolution. U/j — 
and positive initial detuning (ei — S2)/ r f = (6.3 — 6.0)/l = 
0.3-grey dashed line; (7/7 = and negative initial detuning 
{ei~e 2 )h = (6.0-6.3)/l = -0.3-black dashed line; U/j = 3 
and positive initial detuning (ei— £2)/7 = (6.3 — 6.0)/l = 0.3- 
grey line; U = 3 and negative initial detuning (ei — £2)/"/ = 
(6.0 - 6.3)/l = -0.3-black line. 



IV. CHARGE RELAXATION IN THE 
COUPLED QDS IN THE PRESENCE OF 
ON-SITE COULOMB REPULSION 

We start by discussing the case of weak Coulomb in- 
teraction when the ratio Uri2(t)/j < Ae/j is fulfilled 
(FigJSJ). If the initial detuning between energy levels has 
positive value (ei > £2 FigOk) localized charge relaxation 
rate in the first QD and full charge density in the sec- 
ond QD increase (Fig|S^,,b grey line) in comparison with 
the case when Coulomb interaction is absent (Figj5^i,b 
black line). Relaxation rate increases due to the decreas- 
ing of initial detuning value Ae/7 caused by Coulomb 
interaction (FigJSh., Fig0]). Fig 0] demonstrates the de- 
tuning time evolution and reveals that in the absence of 
Coulomb interaction energy levels detuning has constant 
value (FigU]grey dashed line). The presence of on-site 
Coulomb repulsion results in the smaller detuning values 
in the system except the initial time moment and time 
period when the system is quite empty and Coulomb in- 
teraction can be neglected (FigJUgrey line). 

In the opposite case of negative initial energy levels 
detuning (ei < e 2 Fig|3b) even small values of Coulomb 
interaction results in the effective increasing of energy 
levels spacing (Fig [J] black line). Consequently, relax- 
ation rate in the first QD and full charge density in the 
second QD decrease (FigfS^, dashed line) in comparison 



with the case when Coulomb interaction is absent (FigfSjD 
dashed line) and detuning takes on constant value. 

Let us now focus on the charge relaxation processes 
due to the presence of strong Coulomb interaction 
(C/ri2(i)/7 > Ae/7) and positive initial detuning (FigJS]). 
In this case filling numbers time evolution in the first 
QD reveals three typical time intervals with different val- 
ues of the relaxation rates. The first one corresponds to 
the time interval < t < t^max-, where t^max- is an 
instant of time when the filling numbers in the second 
QD reach maximum value n2max (FigHb)- Simultane- 
ously filling numbers in the first QD demonstrate the 
relaxation rate changing (bend formation) (FigJBJi). The 
appropriate detuning behavior is presented on the FigJT] 
One can find that the presence of strong on-site Coulomb 
repulsion results in fast compensation of initial positive 
detuning. Further time evolution leads to the forma- 
tion of negative detuning between energy levels and it 
is clearly evident that the filling numbers maximum in 
the second QD corresponds to the maximum energy lev- 
els spacing. This time interval reveals charge relaxation 
with the typical rate very close to the 7 re s- Coulomb 
interaction increasing in this time interval results in the 
weak decreasing of relaxation rate in comparison with 
the case when Coulomb correlations are neglected. 

The next time interval t 02mal < t < t i min {t 01mm -is 
an instant of time when the filling numbers in the first 
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FIG. 5: Filling numbers evolution in the first a). m(i) and second b). n%(t) QDs. U/'j = 0-black line both for positive 
(ei — £2)/"/ = (6.3 — 6.0)/l = 0.3 and negative (ei — £2)h — (6-0 — 6.3)/l = —0.3 detuning, [7/7 = 3 and (ei — £2)h = 
(6.3 - 6.0)/l = 0.3 -grey line, U/j = 3 and (ei - e 2 )/j = (6.0 - 6.3)/l = -0.3-dashed line. Parameters T/7 = 0.6 and 7 = 1.0 
are the same for all the figures. 



QD n\{t) achieve minimum value ni m i n ) reveals local- 
ized charge relaxation with the typical rate very close 
to the ^nonres- Filling numbers time evolution in the 
first and second QDs simultaneously demonstrates dip's 
formation (Fig|BJi,b) which corresponds to the local mini- 
mum of the effective detuning (Figj7j. This phenomenon 
can be explained by the influence of the following physical 
mechanism: charge redistribution between the QDs by 
means of relaxation channel with the typical relaxation 
rate {\E\ — I) (eg 127)1 in the presence of strong on-site 
Coulomb repulsion. Charge redistribution strongly gov- 
erns the system behavior in the second time interval due 
to the presence of significant real part of the relaxation 
rate i\E\ — E% |) which is proportional to the value Uri2(t) 
(FigJHlgrey line). In the first time interval relaxation oc- 
curs regularly because real part of the relaxation rate 
( I Ei — E% I ) is negligible or nearly absent due to the small 
values of filling numbers in the second QD. Contribution 
from the charge redistribution mechanism becomes cru- 
cial only when filling numbers in the second QD reach 
maximum value. These effect self-consistently influence 
the charge dynamics in the proposed system by means of 
effective detuning changing. 

The third time interval demonstrates charge relaxation 
with the typical rate very close to the 7 res . This inter- 
val exists due to the decreasing of the filling numbers 
amplitude in the second QD. It leads to the increasing of 
the relaxation rate value in comparison with the previous 
time interval. The presence of on-site Coulomb repulsion 
leads to the small detuning values in comparison with the 
second time interval due to the fact that the system is 
quite empty and Coulomb interaction can be neglected 
(Fig J?]). Consequently charge redistribution also can be 
omitted. Therefore with the increasing of time the effec- 
tive detuning aspire to the value without any Coulomb 
correlations in the system. 

Strong Coulomb interaction significantly influence on 
the relaxation processes. To analyze the mechanism of re- 
laxation law modification one have to examine the relax- 



ation exponents evolution, which determine charge relax- 
ation rates changing in each channel of the QDs. More- 
over we shall analyze time evolution of the preexponen- 
tial factors which govern charge re-distribution between 
the relaxation channels. This analysis will be carried out 
for the most interesting case of positive initial detuning, 
when Coulomb interaction leads to the dip's formation. 

Let us start from the analysis of the exponents evolu- 
tion. Their behavior is just the same for the both QDs 
(FigJS]). In the absence of Coulomb interaction second 
channel relaxation rate Ei — E2 always exceeds the first 
channel E\ — E\ relaxation rate. This ratio between re- 
laxation rates also takes place at the initial time moment 
in the presence of Coulomb interaction. Coulomb inter- 
action results in the dip and peak formation in the sec- 
ond and first relaxation channels accordingly (Fig|S^,b). 
First channel relaxation rate maximum value corresponds 
to the second channel relaxation rate minimum value. 
For the large time values evolution laws demonstrate con- 
stant values of relaxation rates for both relaxation chan- 
nels equal to the values obtained without Coulomb inter- 
action. Splitting of the peak in the first relaxation chan- 
nel and dip in the second one can be seen with further 
increasing of Coulomb interaction. Moreover peaks in 
the first relaxation channel correspond to the dips in the 
second relaxation channel and dip in the first relaxation 
channel corresponds to the peak in the second relaxation 
channel (Fig|5b). 

Let us now focus on the preexponential factors (relax- 
ation channel's amplitudes) time evolution in the pres- 
ence of Coulomb interaction. In the second QD preexpo- 
nential factors time evolution is determined by the same 
law (D(t) and F(t) expression 1501) (Fig. [5] grey line). 
Time evolution of the preexponential factors in the first 
QD is quite different (A(t) and C(t) expression [28]) . 

First relaxation channel amplitude always exceeds sec- 
ond relaxation channel amplitude in the first QD and 
both relaxation channels amplitudes in the second QD 
in the absence of Coulomb interaction. This ratio is also 
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FIG. 6: Filling numbers evolution in the first a). m(i) and second b). 712 (i) QDs. (7/7 = 0-grey line, U/y = 6-grey dashed 
line, 17/7 = 12-black dashed line, U/-y = 14-black line. Parameters (e a - e 2 )h = (6.3 - 6.0)/l = 0.3, T/7 = 0.6 and 7 = 1.0 
are the same for all the figures. 
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FIG. 7: Energy levels detuning time evolution in the case of 
positive initial detuning (ei — £2)/7 = (6.3 — 6.0)/l = 0.3. 
U/j — 0-grey line; U/"/ = 6-grey dashed line; U/~y = 12-black 
dashed line; U/j = 14-black. 



valid at the initial time moment in the case of Coulomb 
interaction between localized electrons. 

Relaxation channel's amplitudes time evolution in the 
second QD (D(t) and F(t) in Fig|9}D) demonstrates max- 
imum which corresponds to the minimum for the both 
relaxation channel's amplitudes in the first QD for small 
Coulomb interaction values. With the increasing of 
Coulomb interaction (Fig{9jD) first relaxation channel am- 
plitude in the first QD A(t) and both relaxation channel's 
amplitudes in the second QD (D(t) and F(t)) tend to 
zero, second relaxation channel amplitude in the first QD 
C(t) tends to the unity. Further time evolution demon- 
strates that all channels amplitudes in both QDs turn 
to constant values equal to the values obtained without 
Coulomb interaction. 

Increasing of the Coulomb interaction results in simul- 
taneous peaks formation for all the relaxation channels 
(FigEfc). 



Comparing the obtained results which describe expo- 
nents evolution (charge relaxation rates in each channel) 
and evolution of the preexponential factors (time evolu- 
tion of the each relaxation channel amplitude) one can 
easily reveal that charge redistribution between the re- 
laxation channels in the same QD occurs. At the initial 
moment most part of the charge is concentrated in the 
first relaxation channel of the first QD and later localized 
charge redistributes to the second relaxation channel of 
the first QD. The following time increasing again demon- 
strates charge localization in the first relaxation channel 
of the first QD. Charge in the second QD is equally dis- 
tributed between both relaxation channels. 

Charge relaxation in the presence of Coulomb interac- 
tion in both quantum dots is determined by the charge 
density redistribution among different channels in the 
same QD and by the changing of relaxation rates of 
each channel. So due to Coulomb interaction the leading 
mechanism of non-monotonic charge relaxation in each 
QD is charge redistribution between the channels in a 
separate QD at particular range of the system parame- 
ters. 

The last point of our discussion deals with compar- 
ison between results obtained for suggested model and 
more natural from experimental point of view situation 
when Coulomb interaction is taken into account in both 
QDs. Fig[TU] demonstrates localized charge time evolu- 
tion for the both relaxation channels. Grey line corre- 
sponds to the case when Coulomb interaction exists only 
in the second QD and black-dashed line describes the sit- 
uation when Coulomb interaction is taken into account 
in the both QDs. It is evident that Coulomb interaction 
in the both QDs results in more rapid compensation of 
energy levels detuning for Ae/7 > and doesn't lead 
to qualitative changes of the obtained results. So (as 
it was mentioned above) for simplicity it is sufficient to 
consider Coulomb interaction only in one QD. In any 
case Coulomb interaction effectively changes level spac- 
ing, which control the charge redistribution between the 
relaxation channels in a single QD as well as between two 




FIG. 8: Each channel relaxation rates time evolution for the values of parameters (ej — £2)/! = (6.3 — 6.0)/l = 0.3; T/7 = 0.6; 
7 = 1.0. Ei - -black line; E 2 - TS-dashed black line; Ex - £|-grey line; a). U/-y = 6; b). (7/7 = 12; c). [7/7 = 14. 




FIG. 9: Each relaxation channel amplitude time evolution as a function of time for the values of parameters (ei — £2)/ '7 = 
(6.3 — 6.0)/l = 0.3; T/7 = 0.6; 7 = 1.0. Preexponential factor for the first relaxation channel in the first QD A(t)-b\ack line; 
preexponential factor for the second relaxation channel in the first QD C(i)-dashed black line; preexponential factors for the 
first and second relaxation channels in the second QD D(t) and F(t)-gcey line. a). (7/7 = 6; b). [7/7 = 12; c). U/j = 14. 



QDs. 

A. Conclusion 

We have analyzed time evolution of localized charge in 
the system of two interacting QDs both in the absence 
and in the presence of Coulomb interaction between lo- 
calized electrons within a particular quantum dot. We 
have found that Coulomb interaction modifies the relax- 
ation rates and the character of localized charge time 
evolution. It was shown that several time ranges with 
considerably different relaxation rates arise in the system 
of two coupled QDs. We demonstrated that the presence 
of Coulomb interaction leads to strong charge redistribu- 
tion between different relaxation channels in each QD. So 
we can conclude that non-monotonic behavior of charge 



density is not the result of charge redistribution between 
the QDs but is determined by charge redistribution be- 
tween the relaxation channels in a single QD. 

In any real situation Coulomb interaction is present in 
both QDs. But for simplicity it is sufficient to consider 
Coulomb interaction only in one QD since the main role 
of Coulomb interaction is charge redistribution between 
the relaxation channels in a single QD as well as between 
QDs due to modification of the detuning between energy 
levels in QDs. 
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